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BEAM  TRAPPING  IN  A  MODIFIED  BETATRON  WITH  TORSATRON  WINDINGS 


I.  Introduction 

1-20 

Extensive  theoretical  studies  over  the  last  few  years 

have  shown  that  the  modified  betatron  accelerator  has  the 

potential  to  generate  very  high  current  electron  beams. 

2 1 

Recently,  a  "table  top"  device  at  the  University  of  California, 

Irvine  has  produced  electron  rings  with  about  200A  circulating 

current.  A  larger  device,  currently  under  testing  at  the  Naval 

22 

Research  Laboratory,  has  been  designed  to  produce  multi- 
kiloampere  electron  rings  and  to  address  the  critical  physics 
issues  of  the  modified  betatron  concept. 

A  disadvantage  of  the  modified  betatron  accelerator  is  its 

2 

sensitivity  to  the  energy  mismatch.  Whenever  the  energy  of  the 
electron  beam  is  not  precisely  matched  to  the  betatron  magnetic 
field,  the  center  of  gyration  or  equilibrium  position  of  the  beam 
in  the  transverse  to  the  minor  axis  plane  is  shifted  radially  and 
thus  the  probability  for  the  beam  to  strike  the  wall  is 
increased.  To  alleviate  the  energy  mismatch  problems,  a  modified 
betatron  accelerator  should  be  carefully  designed  to  provide  both 
an  accurate  and  precise  betatron  field  and  the  injector 
accelerator  should  generate  electron  pulses  having  energy  that 
precisely  matches  this  field. ^ 

To  reduce  the  orbit  sensitivity  of  the  modified  betatron  to 

2  3 

the  energy  mismatch,  Roberson  et  al.  have  suggested  the  use  of 

a  strong  focusing  field  generated  by  i  =  2  stellarator 

windings.  In  more  recent  studies,  we  have  demonstrated 

numerically  that  the  energy  bandwidth  of  a  modified  betatron 

configuration  assisted  by  a  strong  focusing  field  generated  by 

torsatron  windings  can  be  very  wide.  In  addition,  we  have  shown 

that  the  torsatron  windings  also  substantially  improve  the 

24 

current  carrying  capabilities  of  the  device  and  could  alleviate 
the  beam  displacement  problem  associated  with  the  diffusion  of 
the  self  magnetic  field. 
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However,  the  usefulness  of  the  twisted  windings  in  improving 
the  performance  of  the  modified  betatron  remained  questionable, 
because  a  scheme  to  trap  the  beam  in  such  a  configuration  was  not 
available.  By  reducing  the  orbit  sensitivity  to  the  energy 
mismatch,  the  strong  focusing  field  made  the  beam  trapping  into 
such  devices  more  difficult. 

2  5 

Recently  Sprangle  and  Kapetanakos  have  developed  a  scheme 
for  trapping  a  beam  into  a  rebatron  accelerator.  This  device  has 
a  magnetic  field  configuration  that  is  similar  to  that  of  a 
modified  betatron  supplemented  with  torsatron  windings.  The 
suggested  scheme  is  based  on  the  dissipative  force  generated  by 
the  resistive  wall  surrounding  the  electron  ring.  Although  such 
a  trapping  scheme  is  appropriate  for  a  rebatron,  it  is  of 
doubtful  usefulness  for  the  modified  betatron.  The  reason  is 
that  in  the  rebatron  the  particle  acceleration  occurs  very 
rapidly,  i.e.,  over  a  few  usee  and  therefore  the  cyclotron 
frequency  corresponding  to  the  betatron  magnetic  field 
(Qzo/y)  very  rapidly  with  time.  Therefore,  the  resistive  wall 
mode  may  not  impose  serious  limitations  on  the  stability  of  the 
ring.  However,  this  is  not  the  case  in  the  modified  betatron, 
because  the  acceleration  occurs  slowly,  i.e.,  over  a  few  msec. 

In  this  paper  we  propose  a  trapping  scheme  that  is 

compatible  with  the  slow  acceleration  of  the  modified  betatron. 

The  proposed  trapping  approach  is  based  on  a  low  amplitude, 

rapidly  varying  magnetic  field  that  shifts  the  beam  equilibrium 

position  near  the  minor  axis  of  the  torus  in  time  shorter  than 

the  beam  bounce  period.  Since  the  new  trapping  scheme  is  similar 

13  14 

to  that  presently  used  in  the  NRL  modified  betatron,  ’  its 
implementation  would  require  only  minor  modifications  in  the 
existing  device. 

Before  addressing  the  beam  trapping,  we  will  derive  two 
equations  that  describe  the  nonlinear  motion  of  the  guiding 
center  of  the  reference  electron  that  is  located  at  the  beam 


centroid.  Using  these  equations  we  will  obtain  a  nonlinear 
expression  that  predicts  very  accurately  the  beam  equilibrium 
position,  even  for  large  displacement  from  the  minor  axis  of  the 
torus.  In  addition,  by  integrating  analytically  the  guiding 
center  equations  we  will  obtain  a  constant  of  the  motion  that 
predicts  accurately  the  topology  of  the  guiding  center  orbits. 
These  analytical  results  have  provided  an  invaluable  insight  in 
the  development  and  understanding  of  the  trapping  scheme. 

II.  Guiding  Center  Motion 

In  this  Section  we  derive  the  equations  describing  the 
guiding  center  motion  for  the  reference  electron  that  is  located 
at  the  centroid  of  a  circular  cross  section  electron  beam.  In 
the  system  of  coordinates  shown  in  Fig.  1,  the  i nstanteneous 
position  of  the  reference  electron  is 

r  =  rQ  +  X  +  ^r  ,  (1) 

z  =  Z  +  5z  ,  (2) 

where  rQ  is  the  major  radius  of  the  torus,  X  and  Z  is  the 
position  of  the  guiding  center  relative  to  the  minor  axis  and 
-5r  and  sz  give  the  position  of  the  gyrating  electron  relative 
to  the  guiding  center. 

Introducing  the  complex  variables 

u  =  r  +  i  z  ,  ( 3 ) 


U 


X  +  i  Z 


(4) 


and  combining  Egs.  (1)  to  (5),  we  obtain 


u=ro=U+6u.  (6) 

Substituting  Eg.  (3)  and  its  time  derivative  into  the  equations 
of  motion 
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where  Q  is  the  cyclotron  frequency  corresponding  to  the 
toroidal  magnetic  field,  v^  is  the  toroidal  velocity,  which  is 
assumed  constant,  E^  and  are  the  image  electric  fields  and 
3^  and  8^  are  the  components  of  applied  and  image  magnetic 
fields.  The  last  terms  consists  of  the  torsatron,  betatron  and 
image  magnetic  fields. 


(a)  Applied  Magnetic  Fields 

The  toroidal  magnetic  field  consists  of  the  applied  toroidal 

field  and  that  generated  by  the  torsatron  windings.  To  lowest 
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order  the  toroidal  field  is  given  by 
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where  B  and  B  are  defined  in  the  local  coordinate  system 

\  A  V  A  J 

ex,  Sy ,  e$  ,  Bq  =  8wI/cL,  I  is  the  current  and  L  is  the 
period  of  the  torsatron  windings, 

0  x  ■* 

B  ?  =  2  n  B  p  K?  (2  a  p  ),  p  is  the  minor  radius  of  the 

S  t  OOc  00 

torsatron  windings,  and  I ^  are  the  modified  Bessel  functions 

and  n  -  2tt/L 

Since  v  is  assumed  constant 
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Ww  *  2ctVq  =  "  2ns/t, 
and  Eq.  (8)  becomes 
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Similarly,  to  the  lowest  order  the  components  of  the 

on 

torsatron  field  transverse  to  the  minor  axis  are 
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Substituting  Eqs.  (11)  and  (12)  into  Eqs.  (13)  and  (14) 


Bt=  8t  cos*  -  8^  sin* 
To  * 


(13) 


Dt  A  at 

B  sin*  +  B  cos*. 


(14) 


and  using  Eq.  (9),  we  obtain 
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In  addition  to  the  torsatron  field,  the  betatron  field 
contributes  to  the  last  term  of  £q.  (7).  Assuming  that  the 
betatron  field  components  vary  as 
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we  obtain 

iBB)  •  13zo  f  1  ♦  „  -  2_  (r  -  1  Z)  1,  (16) 

0 

where  is  the  betatron  field  on  the  minor  axis  and  n  is  the 

external  field  index. 

(b)  Image  Fields  at  the  Ring  Centroid 

An  accurate  self  consistent  determination  of  self  fields  of 
a  high  current  electron  ring  in  a  modified  betatron  configuration 
that  is  supplemented  by  strong  focusing  is  very  difficult, 
because  the  minor  cross  section  of  the  ring  varies  along  the 
toroi dal  di recti  on  . 

However,  here  we  are  interested  in  the  macroscopic  motion  of 
the  ring  and  therefore  on  the  image  fields  that  act  on  the  ring 
centroid.  These  fields  are  not  sensitive  to  the  detailed  shape 


of  the  ring  cross  section  and  thus  can  be  computed  approximately 
assuming  that  the  ring  has  a  circular  cross  section  that  is 
uniformly  filled.  Neglecting  toroidal  corrections,  the  fields  at 
the  reference  particle  are 
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where  r^  is  the  minor  radius  of  the  ring,  a  is  the  minor  radius 
of  the  oerfectly  conducting  torus  and  n  is  the  uniform  electron 
density. 


Since 
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(c)  Taylor  Series  Expansion  of  the  Fields 


To  separate  the  slow  from  the  fast  motion  of  the  reference 
electron  it  is  necessary  to  Taylor  expand  the  fields  about  the 
guiding  center.  Neglecting  terms  higher  than  quadratic,  we  have 
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Eq.  (20)  becomes 
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Furthermore,  we  will  assume  that  the  reference  electron 
rotates  around  the  guiding  center  with  frequency  +  i.e., 


5u  =  Su  t  e  +  ,vj  e 


This  assumption  is  supported  by  the  following  two 

2  4 

observations.  We  have  shown  that  when  a  <<  p  /v,  /■>  is  an 

w  a  w 

ei genf requency  of  the  linear  system  and  that 


We  will  show  later  that  $u  +  ~  I ^  ( 2  a  p)  and 
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in  Eq.  (23)  can  be  neglected,  provided 
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U  <<  v  . 
e 

(d)  Time  Averaged  Equations 


Substituting  into  Eq.  (7)  the  Taylor  expanded  fields  of  Eqs 
(10),  (15),  (16)  and  (19)  using  expression  (21),  omitting  terms 
quadratic  in  3u,  multiplying  the  resulting  equation  by 


e1muwt,  taking  time  average  from  o  to  2n/ow,  and  omitting 


the  small  term  <  U  e1m'Jwt>,  we  obtain  for  m  =  o. 
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In  addition,  for  m=l  the  time  averaged  equation  is 
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Similarly,  for  in  =  -  1 ,  the  time  averaged  equation  gives 
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For  the  parameters  of  interest,  the  underlined  terms  in  Eqs 
(26)  and  (27)  are  small  in  comparison  with  the  non-underl i ned 
terms  and  tnerefore  can  be  omitted. 

In  addition,  assuming  that  y/y  <<  <ow ,  ,  Eqs.  (26)  and  (27) 
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Substituting  Su+  from  (23)  and  (29)  into  the  equation  for  the 
slow  motion  [  Eq.  (25)  ],  we  obtain 
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These  are  the  nonlinear  equations  that  describe  the  guiding 
center  motion  for  the  reference  electron.  The  guiding  center 
gyrates  around  the  equilibrium  position  Xq,  which  can  be 

•  • 

determined  from  (30)  and  (31)  by  taking  X  =  Z  =  o  and  is  given 
by  the  expression 
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(e)  Integration  of  Nonlinear  Equations 


The  nonlinear  equations  (30)  and  (31)  can  be  integrated, 
provided  that  the  spatial  variation  of  n  and  n  can  be 

ys  9  W 

neglected,  i.e.,  when  o  =  qA  =  constant  and  q  -  Q  =  constant 

e  Qo  w  wo 

•  • 

Multiplying  Eq.  (30)  by  X  and  (31)  by  Z  and  taking  their 
difference,  we  obtain 
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The  integration  of  the  last  four  terms  in  Eq.  (33)  is 
elementary.  To  integrate  the  first  term,  we  write 
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where  ^  is  a  function  to  De  determined. 
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the  function  ^  is  obtained  from 


=  riQ  (2<*R) 
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Similarly,  the  second  term  gives 


'l>2  =  rI2  (2ctR)  1  3  (  2«R  )  dR  +  fl  4  (  2aR  )  I  3  (  2aR  )  dR=  l  3  /  Za  (35b) 


Substituting  Eqs.  (34)  and  (35)  into  E  q .  (33)  and  carrying 
out  the  time  integration,  we  obtain 
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where  K  is  a  constant,  0.  =  -  (  q  +  q  )  =  constant  and 
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Equation  (38)  describes  the  nonlinear,  slow  motion  of  the 
beam  centroid  in  the  transverse  to  the  minor  axis  plane,  the 
constant  K  is  determined  from  the  initial  conditions  and  in  order 
the  orbit  to  pass  through  the  minor  axis  (  X  =  Z  =  o)  K  should 
be  zero. 

Near  the  minor  axis,  i.e.,  when  2a  R  <<  1  and  R/  <<  1, 

a 


Eg.  (38)  can  be  linearized  and  the  resulting  expression, 
when  the  small  term  that  contains  I o  is  omitted,  is 


(39; 
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The  displacement  of  the  orbit  center  from  the  minor  axis  :ar 
be  easily  determined  from  (39)  and  is  given  by 
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^zo  ro 

where  yq  =  - -  .  and  5y  =  y-yq  is  the  energy  mismatch. 

so 


Equation  (38)  is  plotted  in  Fig.  2a  for  various  values  of 
constant  K,  for  zero  electron  beam  density.  The  various 
parameters  for  the  results  of  Fig.  2  are  listed  in  Table  I. 

The  number  marked  in  every  fourth  macroscopic  beam  orbit  is 

3 

equal  to  10  K.  Orbits  shown  with  solid  lines  correspond  to  K 
that  is  greater  than  zero  and  those  shown  with  a  dashed  line 
correspond  to  K  that  is  less  than  zero. 

All  the  orbits  close  inside  the  vacuum  chamber.  However, 
those  of  them  that  lie  less  than  a  minor  ring  radius  from  the 
wall  are  not  useful,  because  when  the  beam  is  moving  along  one  of 
these  orbits  it  will  scrape  the  wall. 


Figure  2  b  shows  an  electron  orbit  obtained  from  the 
numerical  integration  of  orbit  equations  superimposed  on  a 
macroscopic  orbit  of  Fig.  2a  that  corresponds  to  the  same 
conditions.  In  both,  the  numerically  obtained  orbit  and  the 
macroscopic  orbit  from  the  constant  of  the  motion  the  toroidal 
corrections  in  the  torsatron  field  have  been  omitted.  It  is 
quite  obvious  that  the  constant  of  the  motion  predicts  accurately 
the  guiding  center  motion.  The  period  of  the  particle  gyrating 
around  the  guiding  center  is  about  1  nsec,  in  agreement  with  Eq. 
(22)  and  the  radius  of  gyration  decreases  as  the  electron  moves 
toward  the  minor  axis,  in  agreement  witn  Eqs.  (28)  and  (29). 

The  predictions  of  the  constant  of  the  motion  remain  rather 
accurate  even  when  the  toroidal  corrections  of  the  torsatron 
field  are  included  in  the  orbit  equations.  Results  are  shown  in 
Fig.  3c . 

For  the  values  of  the  parameters  listed  in  Table  I,  Eq.  (32) 
predicts  that  the  nonlinear  equilibrium  position  is  located  at 
X  ~  10.2  cm.  To  determine  accurately  tne  equilibrium  position, 
we  have  repeated  the  computer  run  of  Fig.  2a  using  a  finer 
scale.  The  results  are  shown  in  Fig.  2d  and  indicate  an 
equilibrium  position  that  is  very  close  to  10.2cm.  To  check 
these  predictions,  we  made  a  computer  run  with  the  electron 
initially  placed  at  X  =  10.0  cm,  Z  =  o.  We  have  observed  that 
the  electron  started  to  gyrate  around  the  X  =  10.2b  cm,  Z  *  o 
point  with  a  radius  of  a  few  millimeters.  This  test  demonstrates 
unequivocally  that  Eq.  (32)  predicts  accurately  the  equilibrium 
position  even  when  the  parameter  2  R  >  1. 

However,  this  is  not  the  case  with  the  linearized  expression 
given  by  Eq.  (40).  For  the  parameters  of  Table  I  Eq.  (40) 
predicts  a  .i r  =  37  cm,  i.e.,  it  overestimates  the  equilibrium 
position  by  about  a  factor  of  four. 
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(f)  Toroidal  Corrections 


So  far,  we  made  the  assumption  that  the  aspect  ratio  a/rQ  is 
very  small  and  thus  all  the  toroidal  corrections  associated  with 
the  fields  have  been  neglected. 

For  devices  with  parmeters  similar  to  those  listed  in  Table 
I  and  II  the  toroidal  corrections  associated  with  the  torsatron 
fields  change  only  slightly  the  ring  orbits  as  it  is  manifested 
by  the  results  of  Fig.  2c.  Even  with  the  torsatron  field 
corrections,  we  were  able  M  derive  the  two  slow  equations  of 
motion,  but  we  were  unable  to  analytically  integrate  them  and 
thus  to  derive  a  constant  of  the  motion. 

In  contrast  to  the  torsatron  field,  the  toroidal  corrections 
of  the  self  fields  become  very  important  when  the  ring  current 
reaches  a  few  kiloamperes.  In  addition  to  changing  the  speed  of 
rotation  around  the  equilibrium  position,  the  self  field  toroidal 
corrections  modify  the  value  of  the  vertical  field  required  to 
confine  the  beam  at  its  equilibrium  position. 


To  lowest  order,  the  self  field  toroidal  corrections  can  be 
included  in  the  slow  equations  of  motion  by  replacing  the  last 
term  in  Eq.  (30)  by 
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Similarily,  the  last  term  in  Eq.  (32)  must  be  replaced  by 
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Equations  (41)  and  (42)  have  been  de  r i ved  under  the  assumption 


that  y/y  <<  u  .  This  assumption  is  valid  because  y  varies  in  the 


time  scale  of  the  slow  beam  motion.  However,  the  y  in  (41)  and 
(42)  is  not  the  same  with  that  at  the  injection  point.  The 
two  y  s  are  related  by 
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in  Eq.  (43)  tne  variation  of  ■<  associated  with  the  toroidal 
electric  field  that  is  generated  by  the  time  changing  self 
magnetic  vector  potential  as  the  beam  centroid  moves  in  the 
transverse  plane  has  been  neglected. 


III.  Beam  Trapping 


The  minimum  requirements  for  trapping  an  electron  beam  in  a 
toroi dal  devi ce  are: 

1.  During  the  first  revolution  around  the  major  axis,  i.e., 

in  a  time  t  <  (2ir/Q2Qy  )  the  beam  should  drift  along  the 

guiding  center  orbit  a  distance  greater  than  rfa  +  r ,  where 

r .  is  the  minor  radius  of  the  beam  and  r .  .  is  the  injector 

b  i  n  j 

radius. 

2.  In  a  bounce  period,  ig  i.e.,  in  the  time  the  beam 
completes  a  revolution  around  the  equilibrium  postion  the  guiding 
center  orbit  should  be  deformed  radially  near  the  injector  at 
least  as  much  as  the  distance  rb  +  rinj  * 

3.  Within  a  time  interval  of  a  few  bounce  periods  the  beam 
should  drift  near  the  minor  axis  of  the  torus,  and 

4.  the  radius  of  gyration  around  the  guiding  center  should 
be  substantially  shorter  than  the  radius  of  the  guiding  center. 

The  guiding  center  equations  derived  in  the  previous  section 
provide  a  valuable  guide  in  the  development  of  a  realistic 
trapping  scheme  for  the  following  reasons.  The  beam  o^ift 
velocity  along  the  guiding  center  orbit  can  be  optimized  using 
Eqs.  (30)  and  (31).  From  the  shape  of  the  guiding  center  orbit 
that  can  be  determined  from  Eq.  (38)  and  the  bounce  frequency,  we 
can  select  the  various  parameters  of  the  trapping  field  that  will 
radially  deform  the  guiding  center  orbit  by  the  desirable  amount 
so  the  beam  wil  not  strike  the  injector  after  a  bounce  period. 

In  addition,  Eq.  (32)  gives  the  value  of  the  betatron  field 
required  to  drift  the  beam  near  the  minor  axis.  Finally,  the 
electron  gyroradius  around  the  guiding  center  can  be  minimized 
using  Eqs.  (28)  and  (29). 


For  the  parameters  of  Table  I,  the  value  of  the  betatron 
field  to  have  the  equilibrium  position  of  the  beam  on  the  minor 
axis  of  the  torus  can  be  computed  from  Eq.  (32)  and  is 
approximately  48  G.  Similarly,  when  the  vertical  field  is  35G 


the  equilibrium  position  Xq  is  located  10.2  cm  from  the  minor 


axis  of  the  torus  (see  Fig.  2d).  Therefore,  by  increasing  the 
betatron  field  as  shown  in  Fig.  3,  from  35  G  to  48  G  over  a  time 
period  that  is  comparable  to  the  beam  bounce  period,  the  beam 


equilibrium  position  will  move  from  Xq  =  10.2  cm 


to 


Xo  =  0 


Fig.  4  shows  the  orbit  of  the  reference  electron 


for  Bzq  =  48  G,  35G  and  also  for  the  combined  field  plotted  in 


Fig.  3,  This  field  consists  of  a  48G  time  in  dependent  field 
and  a  time  dependent  "trapping"  field  that  varies  as 


t  r 
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for  t  >  t 


t  ro  , 


for  t  <  t 
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where  to  is  the  trapping  initiation  time.  The  values  of  the 
various  parameters  of  the  "trapping"  field  for  the  results  shown 
in  Fig.  4  are  listed  in  table  II. 

Over  the  first  20  nsec,  i.e.,  over  a  revolution  around  the 
major  axis  the  beam  drifts  8  cm  away  from  the  injector,  which  is 
more  than  sufficient.  However,  the  electron  drifts  around  the 
new  equilibrium  position  at  X  =  Zq  =  o  very  slowly.  The  total 
time  required  to  complete  the  bounce  in  Fig.  4  is  in  excess  of 
1.3  y s  e  c . 

To  determine  the  sensitivity  of  the  orbit  to  the  initial 
azimuthal  position  of  the  injector,  we  made  a  series  of  runs  with 
the  injector  always  at  the  same  X  and  Z  but  at  different  S 
positions.  Figure  5  shows  the  results  when  the  injector  is  moved 
from  S  =  L/4  to  S  =  o.  As  shown  in  table  II  the  remaining 


parameters  for  this  run  are  identical  to  those  of  Fig.  4.  It 
should  be  noticed  that  the  final  guiding  center  radius,  i.e., 
when  the  electron  gyrates  around  the  XQ  =  ZQ  =  o  equilibrium 
position  is  appreciable  smaller  in  Fig.  5  than  in  Fig.  4.  The 
reason  is  that  in  Fig.  5  when  the  betatron  field  is  35G  the 
constant  of  the  motion  K  is  zero,  as  manifested  by  the  fact  that 
the  orbit  passes  through  the  minor  axis.  However,  the  constant  K 
for  the  35G  orbit  in  Fig.  4  is  positive.  The  trapping  field 
increases  the  value  of  K  by  about  the  same  amount  in  both  runs. 
Thus,  the  final  K  is  greater  for  the  run  of  Fig.  4  than  for  the 
run  of  Fig.  5.  Since  the  guiding  center  radius  near  the  minor 
axis  scales  [see  Eq.  (39)]  as  the  square  root  of  K,  it  is 
reasonable  to  have  a  larger  radius  in  Fig.  4  than  in  Fig.  5. 

There  is  an  additional  difference  between  the  results  of 
these  two  runs,  which  is  not  apparent  in  the  figures.  In  Fig.  5, 
at  the  injection  point  the  torsatron  field  is  directed  downwards 
and  thus  the  force  -  e  v0  B  *  is  radially  outwards.  As  a 
result,  the  initial  motion  of  the  electron  is  outwards.  In 
contrast,  the  -ev^  force  at  the  injection  point  in  Fig.  4 
is  radially  inwards  and  so  is  the  initial  motion  of  the  electron. 

The  results  shown  in  Figures  4  and  5  were  obtained  for  zero 
electron  beam  current.  However,  the  trapping  mechanism  remains 
applicable  even  at  high  beam  current  1^.  Results  are  shown  in 
Fig.  6  for  I  b  =  10  KA.  The  rest  of  the  parameters  for  this  run 
are  listed  in  Table  III.  At  high  beam  current  the  direction  of 
the  slow  motion  is  opposite  to  that  for  Ib  =  0,  because 
the  f  X  £  drift  exceeds  the  combined  curvature  and  gradient 
drifts.  Furthermore,  in  contrast  to  the  low  current  case  the 
equilibrium  postion  moves  toward  the  major  axis  as  the  betatron 
magnetic  field  is  reduced.  This  result  is  in  agreement  with  Eq. 
(40). 

Figure  7a  shows  the  value  of  the  vertical  magnetic  field 
required  to  confine  the  beam  at  its  equilibrium  position 


(equilibrium  magnetic  field)  when  1^  =  0  and  Fig.  7b  shows  the 
corresponding  results  for  I ^  =  10  KA  and  two  different  initial 
injection  positions  X.jnj  *  0  and  12.5  cm.  The  results  of  Fig.  7b 
have  been  obtained  from  Eq.  (32)  with  the  modification  of 
Eq.  (42).  The  various  parameters  for  this  run  ( I b  =  10  KA)  are 
listed  in  Table  IV.  The  equilibrium  magnetic  field  depends  upon 
tne  injection  position  because,  as  Eq.  (43)  indicates,  y  at  Xq 
is  a  function  of  x-jng*  According  to  Fig.  7b,  when  B 1Q  =  135  G 
the  equilibrium  position  XQ  is  at  112.5  cm.  Similarily, 
when  Bzq  =  119  G  the  equilibrium  position  is  at  XQ  =>  107  cm  and 
finally  when  BZQ  =  104  G  the  equilibrium  position  is  located  near 
the  minor  axis. 

In  Figure  6,  the  betatron  and  trapping  field  have  the  same 
polarity  and  thus  until  the  trapping  is  initiated,  i.e.,  for  the 
first  200  nsec  the  total  vertical  field  is  119  G  and  therefore 
the  equilibrium  position  is  at  107  cm.  For  comparison  the  orbit 
of  the  reference  electron  for  a  time  independent  3ZQ  =  119  G  is 
also  shown  in  Fig.  6.  By  decreasing  the  total  vertical  field 
from  119  G  to  104  G  over  a  time  period  that  is  shorter  than  the 
beam  bounce  period,  the  beam  equilibrium  postion  moves  near  the 
minor  axis  and  the  beam  is  trapped. 


IV.  Conclusions 


For  to w <  (flq/  )»  we  have  derived  two  equations  that  describe 
the  nonlinear,  transverse  motion  of  the  guiding  center  of  the 
reference  electron  that  is  located  at  the  beam  centroid.  Using 
these  equations  we  have  obtained  nonlinear  expression  that 
predicts  very  accurately  the  beam  equilibrium  postion,  even  when 
2  a  R  >  1.  Furthermore,  by  integrating  analytically  the  guiding 
center  equations,  we  have  obtained  a  constant  of  the  motion  that 
predicts  accurately  the  guiding  center  orbits  in  the  transverse 
plane. 

With  the  insight  provided  by  these  analytical  results,  we 
were  able  to  develop  a  realistic  trapping  scheme  that  does  not 
require  the  presence  of  a  resistive  wall  and  therefore  is 
compatible  with  the  relative  slow  acceleration  of  the  modified 
betatron . 

The  proposed  trapping  scheme  opens  the  possibility  of  adding 
strong  focusing  in  a  modified  betatron  accelerator.  Such 
focusing  may  have  several  beneficial  effects  on  the  device, 
because  it 

-  will  alleviate  the  energy  mismatch 

-  could  alleviate  the  beam  displacement  associated  with  the 
diffusion  of  self  magnetic  field 

-  will  improve  the  current  carrying  capabilities  of  the  device 
without  requiring  injection  at  higher  energies. 


Table  I 

Parameters  of  the  runs  shown  in  Fig.  2. 


Torus  Major  Radius  r  Q  ( c  m ) 

= 

100 

Winding  Minor  Radius  pQ(cm) 

13 

Toroidal  Chamber  Minor  Radius  a(cm) 

= 

20 

a  =  2  tt / L  (cm  ) 

= 

0.1 

Field  Strength  Factor  e t 

= 

-0.04275 

Winding  Current  I  ( k  A  ) 

S 

25 

l 

= 

2 

Torsatron  Toroidal  Field  Bsex(kG) 

* 

1.0 

Additional  Toroidal  Field  Bsex  ( k G ) 

* 

3.0 

3etatron  Field  8  z  q  (G) 

= 

35 

Ext.  Field  Index  n 

= 

0.5 

Beam  Electron  Density  nQ(cm  ^ ) 

= 

0 

Initial  y 

a 

2.96 

Table  II 


Fi g.  4  Fi g.  5 


Torus  Major  Radius  rQ(cm) 

— 

10G 

100 

Winding  Minor  Radius  p  ( c  m ) 

- 

18 

18 

Toroidal  Chamber  Minor  Radius  a(cm) 

- 

20 

20 

a  =  2  tt/L  (cm  ) 

= 

0.1 

0.1 

Field  Strength  Factor  e t 

= 

-0.04275 

-0.04275 

Winding  Current  I  ( k  A  ) 

= 

25 

25 

Z 

= 

2 

2 

Torsatron  Toroidal  Field  B$ 

(kG) 

= 

1.0 

1.0 

e-Fold  time  x(nec) 

= 

100 

100 

Trapping  Initiated  at  Time  (nsec) 

= 

100 

100 

Trapping  Field  Amplitude  Btf, 

(G) 

= 

-13 

-13 

Additional  Toroidal  Field  8 

65  (kG) 

= 

3.0 

3.0 

Betatron  Field  Bzq (G ) 

= 

48 

48 

Ext.  Field  Index  n 

s 

0.5 

0.5 

Electron  Beam  Current  I b  (kA) 

= 

0 

0 

Initial  y 

= 

2.96 

2.96 

Initial  Positions  ( cm ) 

X(t=0) 

= 

14 

14 

Z(t=Q) 

= 

0 

0 

S(t=0) 

= 

15.7 

0 

Initial  Velocities 

X ( t  =0  ) 

= 

0 

0 

2  ( t = o ) 

= 

0 

0 

% 

= 

0.941c 

0.941c 

vj 

ii 


l 


I 


29 


Table  III 

Parameters  of  the  run  shown  in  Fig.  6 


Torus  Major  Radius  r0  (cm) 

— 

100 

Winding  Minor  Radius  o  (cm) 

- 

18 

Toroidal  Chamber  Minor  Radius  a  (cm) 

a  =  2  rr/L  (cm  ) 

= 

16 

= 

0.1 

Field  Strength  Factor  et 

= 

-0.04275 

Wi ndi ng  Current  I  (kA) 

= 

25 

9. 

= 

2 

Torsatron  Toroidal  Field  Bs  (kG) 

= 

1.0 

e-Fold  time  t  (nsec) 

= 

100 

Trapping  Initiated  at  Time  (nsec) 

= 

200 

Trapping  Field  Amplitude  (G) 

= 

15 

Additional  Toroidal  Field  B^x  (kG) 

= 

3.0 

Betatron  Field  6 2 Q  (G) 

* 

104 

equilibrium  Mag.  Field  at  Inj.  8  g  ^  (G) 

s 

135  ' 

Ext.  Field  Index  n 

= 

0.5 

Electron  Beam  Current  I b  (kA) 

= 

10 

Initial  y 

= 

2.96 

Minor  Beam  Radius  r^  (cm) 

= 

1.0 

v/y 

= 

0.21 

Initial  Positions  (cm) 

O 

11 

4-> 

X 

- 

13 

Z(t  =  0) 

= 

0 

S(t  -  0) 

= 

15.7 

Initial  Velocities 

X(t  =  0) 

3 

0 

l(t  =  0) 

3 

0 

V, 

= 

0.941c 

v 


Table  IV 

Parmeters  of  the  run  shown  in  Fig.  7 


Torus  Major  Radius  rQ  (cm) 

Winding  Minor  Radius  p  (cm) 

Toroidal  Chamber  Minor  Radius  a  (cm) 
a  =  2  ir/L  (cm  ^ ) 

Field  Strength  Factor 
Wi ndi ng  Current  I  (kA ) 


-0.04275 


Torsatron  Toroidal  Field  Bs  (kG) 
Additional  Toroidal  Field  B®x  (kG) 
Ext.  Field  Index  n 
Electron  Beam  Current  I b  (kA) 
Initial  y 

Minor  Beam  Radius  r  b  (cm) 
v/ y 

Initial  Positions  (cm) 


2.96 


0.21 


Initial  Velocities 


0.941c 


— -160 


Equilibrium 
Position 
X0  =  10.2  cm 


r0  +  X  (cm) 

2(d)  —  As  Fig.  2(b),  but  a  finer  scale 


'Btro  e_(t"to>/T ,  for  t>f 
.Btro  -  for  t<f0 


Betatron,  trapping  and  total  vertical  Held  as  a  function  of  time 


r0  +  X  (cm) 

Fig.  4  —  Orbit  of  reference  electron  at  the  beam  centroid  injected  at  114  cm  from  the  major  axis  when 
B.„  »  35  G  =  constant,  B.0  *  48  G  —  constant  and  in  the  presence  of  the  trapping  magnetic  field. 


The  various  parameters  for  this  run  are  listed  in  Table  II. 


r0  +X0(cm) 

Fig.  7  —  Vertical  magnetic  field  required  to  confine  the  beam  at  its  equilibrium  position  as  a  function 
of  the  equilibrium  position  for  I*  =  0  (Fig.  7a)  and  I4  =  10  KA  (Fig.  7b).  The  various  parameters  for 
these  figures  are  listed  in  Table  IV. 
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Fig.  7  (Corn'd)  —  Vertical  magnetic  field  required  to  confine  the  beam  at  its  equilibrium  position  as  a 
function  of  the  equilibrium  position  for  IA  —  0  (Fig.  7a)  and  I6  =  10  KA  (Fig.  7b).  The  various 
parameters  for  these  figures  are  listed  in  Table  IV. 
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